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I.  .INTRODUCTION 


1 . 1  ffb.iect  and  Scope 

thars  ■  lnve5lrf:gaMeft=-4c~'^e  development  of  a  systematic 
numerical  procedure  for  determining  the  displacements^  strains,  and  stresses 
within  a  plane  continuumy^wnerein  certain  regions  have  heen  strained  heyond  an 
elastic  yield  limit.  Such  a  procedure  should  make  possible  the  observation 
of  the  development  of  the  stress  and  strain  patterns  around  regions  of  high 
stress  intensity,  such_as_j::egioHa— ar.©u.ndr.^eitGh.6B.,»dioles-,--aad.»pG^te4s».©f»usoii«to,, 


The  procedure  is  restricted  to  plane,  static  problems;  and  the 
example  problems  are  further  restricted  to  plane  strain  conditions.  The  pro¬ 
cedure  itself  is  applicable  to  plane  stress  problems  if  the  relations  between 
stress  and  strain  for  plane  stress  conditions  are  substituted  for  those  of 
plane  strain.  The  material  of  the  continuum  is  considered  to  be  isotropic 
and  elastic-perfectly -plastic  and  the  problems  are  solved  for  continuously 
increasing  external  loads.  Unloading  from  a  plastically  strained  state  is 
not  considered. 

The  numerical  procedure  is  essentially  a  relaxation  techniq.ue 
I  ^ 

applied  to  a  discrete  physical  model  composed  of  suitably  arranged  stress 
points  and  mass  points.  Plastic  yie.lding  and  flow  in  the  solid  is  charac¬ 
terized  by  the  corresponding  yielding  and  flow  of  the  stress  points  of  the 
model . ^^ntroduction  of  the  discrete  model  reduces  the  problem  of  the 
continuum  with  an  infinite  number  of  degrees'  of  freedom  to  a  problem  in 
particle  mechanics  with  a  finite  number  of  degrees  of  freedom.  The  primary 
advantage  of  such  a  technique  is  that  it  makes  possible  the  solution  of 
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problems  not  tractable  by  ordinary  mathematical  analysis,  particularly 
problems  involving  partial  loadings  and  complicated  boundary  conditions.  The 
basic  disadvantage  of  the  discrete  model  approach  is  its  very  finiteness — ^ 
stresses  and  displacements  are  defined  only  at  a  finite  number  of  points . 

Hence  frequently  the  finite  model  can  furnish  only  a  rough .quantitative 
measure  of  the  true  but  unknown  solution  in  the  continuum.  To  gain  some  notion 
of  the  accuracy  of  the  model  used  in  this  investigation,  a  problem  in  plane 
elasticity  for  which  there  is  an  analytic  solution  is  solved  by  using  the 
model  and  the  results  compared  to  the  analytic  solution. 

Once  the  level  of  external  loading  has  been  raised  to  a  sufficiently 
high  level,  the  more  highly  stressed  regions  of  the  continuum  begin  to  yield, 
or  flow  plastically.  The  initiation  of  yielding  is  determined  by  the  Misep- 
Hencky  yield  criterion.  Thereafter,  yielded  regions  are  assumed  to  obey  the 
plastic  stress-strain  relations  postulated  by  the  Prandtl-Reuss  theory.  Two 
examples  of  problems  wherein  plastic  flow  has  taken  place  over  a  finite  region 
are  included  to  demonstrate  the  application  of  the  numerical  procedure. 

The  entire  procedure  for  handling  plane  problems  of  contained 
plastic  flow  in  elastic-perfectly -plastic  continua  has  been  coded  for  use  on 
the  IBM  7090  digital  computer.  Only  the  two  numerical  solutions  mentioned 
above  are  included  in  the  thesis;  an  extensive  investigation  of  the  various 
problems  of  Interest  in  contained  plastic  flow  falls  outside  the  scope  of 
this  work. 


1.2  Historical  Notes 

This  brief  review  of  the  literature  is  by  no  means  intended  to  be 
complete.  Only  a  few  of  the  more  important  publications  related  to  the  present 
study  are  discussed.  -Several  of  these  references  (6),  (8),  (15),  (22)^  contain 
1  Numbers  in  parentheses  refer  to  corresponding  entries  in  the  Bibliography. 


extensive  'bi'bliographies  or  footnotes  through  which  more  detailed  access  to 
the  literature  may  be  obtained. 

The  idea  of  replacing  plane  elastic  continua  with  discrete  models 
began  to  attract  researchers'  interest  in  the  early. 1940' s --about  the  same 
time  that  Southwell  (19)  (20)  developed  efficient  and  practical  relaxational 
techniques  for  the  solution  of  highly  complex  systems.  It  is  not  at  all  sur¬ 
prising  that  the  development  of  finite  models  should  have  awaited  more 
efficient  methods  of  computation,  since  by  their  very  nature  solutions  de¬ 
termined  with  the  use  of  models  involve  systems  with  a  large  nimiber  of 
simultaneous  equations.  Hrennikoff  (lO)  and  McHenry  (15)  were  perhaps  among 
the  earliest  of  those  who  Introduced  "frameworks"  or  "lattice  analogies"  to 
solve  problems  in  plane  elasticity.  Using  hexagonal  and  square  patterns  as 
the  basic  module  in  the  discrete  model,  Austin  (3)  and  Dauphin  (5)  made 
informative  comparisons  of  the  model  solution  to  the  exact  analytic  solution 
for  several  problems  in  plane  elasticity.  Wewmark  (15)  gives,  a  good  di'  's- 
sion  of  the  use  of  models  in  several  areas  of  structural  analysis. 

More  recently,  there  has  been  a  renewed  interest  in  the  development 
of  models;  this  is  partly  prompted  by  more  efficient  computational  devices. 
The  advent  of  high-speed  digital  computers  has  induced  many  analysts  to  seek 
discrete  models  suitable  for  digital  computation.  The  work  of  Clough  (4-) 
and  Gaus  (7)  is  typical  of  the  model  approach  now  being  adopted  in  order  to 
so.lve  continuum  problems  on  computers .  It  is  interesting  to  note  that  none 
of  the  writers  above  make  any  mention  of  attempts  to  extend  their  models 
beyond  the  elastic  range.  Schnobrich  (18)  has  indicated  that  considerations 
for  future  extension  into  plastic  behavior  influenced  the  selection  of  his 
model,  though  his  work  presents  only  ela.stic  results. 

The  scientifip  study  of  the  theory  of  plasticity  seems  to  be  much 
older  than  any  serious  study  of  finite  models,  for  it  extends  back  at  least 
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to  CoulomlD  and  his  study  of  yielding  in  soils  in  I773.  Any  number' of  readable 
texts  in  the  elementary  theory.  (8),  (9),  (.I6),  (17)  are  available,  though  the 
presentation  here  follows  most  closely  that  in  Prager  and  Hodge  (17)  and 
Hoffman  and  Sachs  (9).  The  only  successful  numerical  solutions  of  problems 
in  contained  plastic  flow  known  to  the  author  are  those  presented  by,  Allen 
and  Southwell  (l)  and  Jacobs  (ll) .  Their  solutions  are  obtained  by  a  rather 
tedious  manual  relaxation  technique  which  yields  values  of  the  stress  function 
from  which  the  stresses  are  computed. 

In  suimmary.  then,  it  appears  that  both  the  theory  of  plasticity  and 
the  theory. of  models  have  attracted  the  efforts  of  able  researchers,  though 
there  have  been  few,  if  any,  attempts, to  apply  the  theoiy  directly  to  a 
discrete  model.  Accordingly,  it  is  the  purpose  of  this  investigation  to 
develop  a  numerical  procedure  for  solving  problems  of  contained  plastic  flow 
with  the  use  of  a  discrete  model. 

1.3  Notation 

The  following  notation  has  been  adopted  for  use  in  this  thesis. 

X  direction  of  axis 

y  direction  of  axis 

z  direction  of  axis  (perpendicular  to  the  plane  of  analysis) 

u  disp.lacement  in  x  direction 

V  displacement  in  y  direction 

T],h  displacement  in  horizontal  direction 

|,v  displacement  in  vertical  direction 

E  Young's  modulus 

V  Poisson's  ratio 

E 

G  shear  modulus  =  o'h  .,,T 
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K 


B 

E® 


E 

bulk  modulus  = 

yield  stress  in  simple  shear 
total  stress  tensor 
spherical  stress  tensor 
deviator  stress  tensor 
total  strain  tensor 
spherical  strain  tensor 
deviator  strain,  tensor 


mean 


normal  stress  -  ^  (a  +  a  +  o  ) 
5  X  y  z' 


normal  component  of  S  in  x  direction  =  -  s 

normal  component  of  in  y.  direction  =  a  -  s 

y 

normal  component  of  in  z.  direction  =  -  s 

mean  normal  strain  =  ^  (e^  + 
normal  component  of  in  x  direction  =  "  s 

normal  component  of  in.  y,  direction  =  e  -  e 

y 

normal  component  of  in  z  direction  =  -  e 

principal  nqrmal  component  of  stress  deviation  =  -  s 

principal  normal  component  of  stress  deviation. =  Ug  -  s 

principal  normal  component  of  stress,  deviation  =  a_  -  s 

principal  normal  component  of  strain  deviation  =  e 

principal  normal  component  of  strain  deviation  =  -  e 

principal  normal  component  of  strain  deviation  =  -e 

first  invariant  =  s^  +  Sg  +  s^ 

X  2  2  2 

second  invariant  =  t-  (s,,  +  +  s_) 

4  ±  d  ^ 

1  ^  "Z  "Z 

third  invariant  =  —  (s^  +  Sg  +  s^) 


W 


work  performed  by  stresses  during  a  plastic  distortion 
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F 

S 

X 

Y 

I 

A 


xy 


xz 


yz 


axial  force  component  at  a  stress  point 
shear  force  component  at  a  stress  point 
tody .  force  per  unit  voliune 
body  force  per  unit  volume 

moment  of  inertia  of  a  unit  width  of  the  reinforcing  frame 

cross-sectional  area  of  a  unit  width  of  the  reinforcing  frame 

element  of  the  total  stress  tensor 

element  of  the  total  stress  tensor 

element  of  the  total  stress  tensor 

element  of  the  total  stress  tensor 

element  of  the  total  stress  tensor 

element  of  the  total  stress  tensor 

principal  normal  stress 

principal  normal  stress 

principal  normal  stress 

element  of  the  total  strain  tensor 

element  of  the  total  strain  tensor 

element  of  the  total  strain  tensor 


7  element  of  the  total  strain  tensor 
xy 

7  element  of  the  total  strain  tensor 
xz 


7  element  of  the  total  strain  tensor 
yz 


principal  normal  strain 
Sg  principal  normal  strain 
principal  normal  strain 

X.  factor  of  proportionality  between  stress  and  strain  rate, 

horizontal  or  vertical  distance  between  mass  points 
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6  distance  along  x  or  y.axis  iDetween  mass  points 

0^^  level  of  external  load,  at  which  first  yielding  begins 

f  flexibility  coefficient 

P  concentrated  external  load 


I'l.  DESCRIPTION  OF  THE.  MODEL 


2.1  Criteria  for  Selection  .of  a  Model 

Historically,,  there  have  heen  at  least  two  distinct  criteria  for 
selecting  a  finite  mechanical  model  to  replace  a  continue.  Hrennikoff  (lO) 
and  Clough  (k)  hoth  demand  equality  of  deformaitlon  between  model  and  continuimi 
under  similar  loading  conditions..  It  is  interesting  in  this  regard  to  quote 
Hrennikoff  (lO). 

It  is  now  possible  to  formulate  the  basic  principle 
governing  determination  of  the  framework  pattern.  The 
necessary  and  sufficient  condition  for  equivalence  of 
infinitesimal  framework  and  solid  material, is  equality 
in  deformability  of  the  two... 

Hpennlkoff ,!  s  lapplication  of  .this  cfiterioh '  is-  qUest-lOhablb,  .  sihce  sevefdi-.of  his 
simple  framework  patterns,  deform  as  doe.s  the  continuum  only  if  Poisson's  ratio 
has  the  value  l/3.  Michell  (li^)  shows,  however,  that  at  least  for  simply 
connected  regions  the  values  of  the  elastic  constants  do  not  affect  the  compu¬ 
tation  of  the  stresses  if  the  boundary  conditions  are  specified  by  loading 
conditions  rather  .than  by  displacement  conditions.  Nevertheless,  any  criterion 
which  restricts  the  value  of  a  material  constant  to  a  specific  value  cannot 
be  completely  satisfactory .for  treatment  of  the  most  general  problems. 

A  second  criterion  that  is  sometimes  used  in  the  selection,  of  a 
model  was  mentioned  by  Newmark  (I5)  and  attempted  by  Gaus  (7),  and  was 
explicit.ly  proposed  by.Ang  (2)  in  the  development  of  the  model  which  is  used 
in  this  thesis.  The  criterion  is  that  there  be  a  mathematical  consistency 
between  the  finite  equations  governing  the  behavior  of  the  model  and  the  dif¬ 
ferential  equations  governing  the  behavior  of  the  continuum.  By  this  is  meant 
that  the  equations  for  strains,  stresses,  equilibrium,  and  compatibility  which 
are  derived  directly  from  the  model  should  be  the  same  as  a  set  of  finite 

difference  equations  of  the  corresponding  differential  relations,  governing  the 
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continuum.  If  this  requirement  is  met  the  requirement  of  equal  deformahility 
of  a  model  and.  the  corresponding  continuimi  will  automatically  he  satisfied,  . 
and  no  restriction  need  he  placed  on  the  value  of  Poisson's  ratio  or  of  any 
other  elastic  constant. 

2.2  Description,  of  the  Model 

The  model  used  in  this  investigation  possesses  all  the  requirements 
of  the  second  criterion  cited  above.  The  essential  characteristics  of  the 
model  are  shown  in  Fig.  1,  wherein  a  square  grid  has  been  superposed  on  the 
continuum.  The  mass  of  the  continuum  is  concentrated  at  the  intersections 
of  the  grid  lines.  Each  of  the  mass  points  is  connected  through  stress  points 
to  the  neighboring  mass  points.  Three  components  of  stresses  and  strains  are 
defined  at  each  stress  point  (two  perpendicular  axial  components  and  a  shear 
component) .  Displacements  in  the  continuum  are  defined  only  at  the  mass 
points  while  stresses  and  strains  are  defined  only  -at  the  stress  points. 
Modifications  of  the  model  to  include  a  stiffened  rectangular  opening  are  also 
shown  in  Fig.  1. 

There  are  two  important  advantages  of  the  model  configuration 
described  above.  First,  all  elements  of  the  strain  tensor  and  the  stress 
tensor  are  defined  at  the  same  point.  This  is  an  important  characteristic 
of  the  model,  especially  in  extending  its  use  to  problems  of  plasticity. 
Second,  horizontal  and  vertical  boundaries  of  the  model  contain  only  mass 
points.  Thus  boundary  conditions  given  in  terms  of  either  external  tractions 
acting  on  the  mass  points  or  prescribed  displacements  of  these  mass  points 
can  be  applied  vrith  equal  ease. 
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2 . 3  Relation  of  the  Model  to  Finite  Difference  Equations  of  the ■ Continuum 

The  material  presented  in  this  section  follows  closely  that  given 
by  Ang  (2).  For  purposes  of  illustration  the  following  notation  is  used. 
Superscript  letters  refer  to  stress  point  locations.  Subscript  letters  x 
and  y  refer  to  the  directions  of  the  axes.  Subscript  numbers  refer  to  mass 
point  locations.  Displacement  components  in  the  x  and  y, directions  are  given 
by  u  and  v,  respectively.  Sign  convention  is  that  established  by  Timoshenko 
(21)., 

The  components  of  the  strains  at  a  typical  stress  point  "a"  are 
defined^  with  reference  to  Fig.  1^  as  follows; 


Ur-,  -u,  , 
a  _  3k  k3 

X  & 

a  _ 

V  8 

a  _ 

'^xy  6  6 


(1) 


These  strains,  which  are  derived  directly  from  the  model,  are  identical  to 
the  finite  difference  expressions  for  the  differential  strain-displacement 
relations  of  the  classical  theory  for  plane  continue  under  small  deformations 


5u 

^x  ~  dx 
&v 

“  dy 

y  _  Su  ^  5v 
'^xy  5y  Sx 


(2) 


The  equation  of  equilibrium,  in  the  x  direction,  for  a  typical 
Interior  mass  point  at  ■"l*-3"  is,  (see.  Fig.  2) 
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where  X  is  the  "body  force  per  unit  volume,  The  yoluine  of  a  parallelepiped 

2  8^ 

of  unit  thickness  and  area  X  f=  is  considered  concentrated  at  each  mass 
point,  If  the  thickness  of  the  model  is  taken  as  unity  in  the  z  direction, 
forces  at  the  stress  point  "a"  are  obtained  from  the  stresses  as  follows : 


Using  .Eqs,  ()+)  in  Eq,  (j),  the  following  equation  of  equilihrium^  in  terms 
of  stresses ;  is  obtained; 


a 

a 

X 


c 

a 

X 


8 


+ 


X  =  0 


(5a) 


A  similar  equation  is  obtained  in  the  y  direction: 


+  Y  =  0 


(513) 


These  equilibrium  equations,  (5a)  and  (5b),  are  identical  to  the 
finite  difference  expressions  for  the  differential  equations  of  equilibrium 
governing  the  corresponding  continuum; 


St 

+  X  =  0 

oy 


9a  9i: 

oy  ox  • 


(6a) 


(6b) 
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The  strains  in  the  model  will  necessarily  satisfy  the  compatibility 
relation,  since  strain  compatibility  is  essentially  a  requirement  placed  pn 
the,  three  components  of  strain  in  order  to  insure  that  the  three  strain 
components  correspond  to  a  physically- possible  displacement  configuration, 
fhe  model  deals  directly  with  displacements,  and  the  strains  are  defined 
directly  in  terms  of  these  displacements.  Hence,  it  can  be  expected  that  the 
strains  derived  from  the  displacements  of  the  model  will  identically  satisfy 
the  compatibility  condition? 

It  is  also  possible  to  express  the  equations  of  equilibrium  in  terms 
of  displacements.  This  is  done  below  for  a  linearly  elastic  solid  in  plane 
strain.  Similar  relations  exist  for  plane  stress  conditions.  For  this 
purpose  it  is  necessary  to  express  the  three  force  components  at  the  stress 
point  "a"  in  terms  of  displacements,  as  follows: 


_a  _  B 

X  (1-i-v)  (1^2v) 


(1.,) 


„a  _  .  E.  . 

^y  -  (l:fy)(l-2v) 


53' „  ^5h-^h3 
6 


•(i-v)-  ■  +  V 


8  ■ 
2 


xy  2(l-t-v) 


8  8 


(T) 


Eqs..  (7)  are  essentially  Hooke's  stress-strain  relationships  for  plane  strain. 
Substitution  of  these  and  similar  relations  for  the  forces  originating  at  the 
other  stress  points  intoEq.  (5)  results  in  the  following  equation  of  equi¬ 
librium  in  the  x  direction,  in  terms  pf  displacements: 


E 


(l-fv)(l-2v) 


u^,  -2u,  ,+u,„ 

2(l-v) 

8 


(l.2v)  ^ 

8^  8^ 


rt-  X  F  0 


(8) 
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A  similar  equation  exists  for  equilibrium  in  the  y  direction.  Note  that  this 
equation  is  identically  the  same  as  a  finite  difference  equation  for  the 
differential  equation  of  equilibrium,  .Eq.  (.9) ,  governing  the  continuum. 


(i+v)(1^2v) 


2(l-v) 


3x2 


(l-2v) 


S^v-^ 

3x3y. 


+  X  =  0 


(9) 


2,h-  Boundary  Conditions 

In  general,  boundary  conditions  (for  either  continue  or  discrete 
models)  can  be  of  two  types:  either  the  forces  acting  along  some  boundary 
or  the  displacements  on  the  boundary  are  prescribed.  As  pointed  out  eaplier, 
either  type  of  condition  can  be  imposed  on  the  model.  A  few  examples  are 
given  below  to  indicate  how  boundary  conditions  are  prescribed  for  the  model. 

■  For  greater  flexibility  and  ease  in  programming,  an  extra  line  of ' 
mass  points  has  been  included  on  each  of  the  four  sides  of  the  rectangular 
model,  as  indicated  in  Fig.  1  by  dotted  lines.  Thus  if  the  continuum  being 
simulated  is  to  be  ten  K  high  and  eight  X  wide  (demanding  a  grid  of  eleven 
rows  and  nine  columns),  there  will  actually  be  thirteen  rows  and  eleven 
columns  in  the  complete  description  of  the  model.  Suppose  that  the  continuum 
is  known  to  possess  symmetry  about  a  vertical  axis  through  a  column  of  mass 
points.  The  boundary  condition  on  the  right  edge  of  the  model  (see  Fig.  l) 
is  specified  as  follows: 

"^16  =  ""il^ 

i  =  1,  2,  ...,  6  (10) 

If  an  external  load  is  to  be  applied  to  the  top  surface  of  the 
continuum,  the  model  will  have  equivalent  concentrated  loads  applied  at  the 
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second  row  of  mass  points,  and  the  extra  top  row  of  mass  points  will  he 
neglected  completely.  If  it  is  desired  to  hold  the  hase  of  the  continuum 
fixed  against  displacement, -'.the  bottom,  row  of  extra  mass  points  is  simply 
given  a  zero  displacement. 

■In  problems  for  which  the  model  is  being  used  to  simulate  an 
infinite  half-space,  the  problem  of  what  boundary  conditions  to  impose  on, 
the  left-most  column  of  extra  mass  points  arises.  For  vertical  loadings 
which  are  symmetric  about  the  center  line,  it  has  been  assumed  that  the 
horizontal  displacements  of  this  left-most  column  are  zero  and  that  the 
vertical  displacements  of  this  left-most  column  of  mass  points  will  be  equal 
to  the  vertical  displacements  of  the  column  of  mass  points  immediately  to 
the  right  of  this  boundary  column.  When  these  vertical  and  horizontal  motions 
are  resolved  into  displacements  in  the  x  and  y  directions,  the  boundary  con¬ 
ditions  become 


u.T 

il 


V 


il 


2  ^^^12  +  ^2) 


2  ^^12  ^  ^i2> 


i  =  1,  .  2,  . . . ,  6 


(11) 


These  examples  indicate  the  manner  in  which  boimdary  conditions 
are  prescribed  for  the  model.  A  variety  of  practically  significant  condi¬ 
tions  can  be  conceived,  and  several  different  sets  of  boundary  conditions 
were  actually  Investigated  during  preparation  of  the  numerical  examples .  An 
extensive  treatment  of  the  effect  of  various  boundaiy  conditions  on  the  stress 
and  displacement  patterns  is  beyond  the  scope  of  the  present  work. 


2.5  ■  Modification. of  the  Model  to  Include  Interior  Rectangula,r  Openings 

An  example  of  the  adaptability  of  the  model  approach  to  structural 
analysis  is  given  in  the  problem  of  determining,  the  stress  pattern  within  a 
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plane  solid  around  a  rectangular  opening,  which  opening  may  or  may  not  he 
reinforced.  In  the  general  case,  it  is  supposed,  that  the  opening  is 
reinforced.  If  the  opening  is  to  he  a  cavity, only,  the  modulus  of  elasticity 
of  the  reinforcing  material  is  set  equal  to  zero. 

The  reinforcement,  if  any,  in  the  continuum  is  replaced  in  the  model 
hy  a  system  of  moment  and  axial  springs .  As  shown  in  Pig,.  1,,  the  moment 
springs  are  located  at  the  mass  points,  and  the  axial  springs  span  from  mass 
point  to  mass  point  in  either  a  vertical  or  a  horizontal  direction.  By  means 
of  the  moment  springs,  shear  forces  due  to  moments  in  the  reinforcing  continutmi 
can  he  simulatedj  axial  springs  simulate  the  direct  tensile  or  compressive 
forces  in  the  reinforcing  continuum.  Tensile  forces  in  the  axial  springs  are 
taken' as  positive.  Sign  convention  for  positive  moments  is  basically  dictated  . 
hy  the  requirements  for  positive  shears  arising  as  a  result  of  these  moments. 
This  sign  convention  is  shown  in  Pig.  J. 

The  vertical  or  horizontal  shearing  forces  acting  on  each  mass 
point  (depending  on  whether  the  mass  point  is  on  a  horiziontal  or  vertical 
reinforcing  section,  respectively)  can  he  calculated  from  the  differences  in 
the  moments  acting  at  three  consecutive  mass  points.  Until  a  moment  spring 
begins  to  yield,  the  moment  can  he  computed  directly  from  the  displacements 
(Pig.  1)  as  follows: 


M. 


-  El 


2(Uj5+Vj^)  +  (Ujg+Vjg) 


(12) 


where  E,  is  the  modulus  of  elasticity  of  the  reinforcing  material,  modified 
for  plane  strain,  and  I  is  the  moment  of  inertia  of  a  unit  width  of  reinforce¬ 
ment.  After  a  moment  spring  has  reached  its  yield  limit,  it  is  ass\mied  to 
hold  the  yield  moment,  even  though  the  rotation  of  the  section  may  increase 


cons Iderahly . 
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The  vertical  or  horizontal  axial  forces  acting  on  each  mass  can  he 
determined  as  the  algebraic  difference  of  the  axial  springs  acting  on  each 
side  of  the  mass  point.  Until  an  axial  spring  yields,  the  force  in  a  single 
axial  spring  can  he  computed  from  the  displacements  (Fig.  as  follows; 


F 


15 


AE 

s/2X 


{U36-V36)  -  ('"35'^35^ 


(13) 


where  A  is  the  cross-sectional  area  of  the  reinforcement,  E  is  the  modulus 
of  elasticity  of  the  reinforcing  material,  modified  for  plane  strain  con¬ 
ditions.  After  an  axial  spring  has  reached  its  yield  limit,  it  is  assumed 
that  the  axial  force  maintains  this  yield  level  regardless  of  the  values  of 
the  surrounding  displacements. 

Once  the  horizontal  and  vertical  forces  acting  on  a  mass  point  as 
a  result  of  the  reinforcement  are  determined,  they  are  resolved  into  x  and 
y  components  and  handled  in  the  same  way  as  the  forces  in  the  rest  of  the 
solid. 

It  is  evident  that  a  mass  point  which  lies  on  an  interior  opening 
will  have  forces  acting  on  it  that  are  different  from  the  forces  acting  on 
a  general  interior  mass  point.  It  is  also  evident  that  the  forces  acting  on 
mass  points  which  lie  on  the  opening  will  vary  depending  on  whether  the  mass 
point  is  on  the  top,  bottom,  side,  or  corner  of  the  opening.  For  this  reason 
a  set  of  operators  has  been  developed  which  computes  the  forces  acting  on  a 
mass  point,  given  the  location  of  the  mass  point. 


III.  CONSTITUTIVE  BEIATIONS  PROM  THE  THEORY  OF  PERFECT  PLASTICITY 


3.1  General  Remarks  on  the  Theory  and  Its  Limitations 

Any  constitutive  relationships  of  the  theory  of  plasticity  may  he 
divided  into  the  following  three  parts; 

(1)  stress-strain  relations  for  the  elastic  region, 

(2)  yield  criterion  to  define  the  initiation  of  yielding, 

(3)  stress-strain  relations  for  the  plastic  region. 

These  three  major  divisions  of  the  theorjr  will  be  discussed  in  turn,  after 
the  associated  assumptions  and  limitations  are  listed  and  after  a  set  of 
notation  that  will  be  useful  in  the  discussion  of  the  theory  is  introduced. 

There  are  three  main  assumptions  underlying  the  theory  of  perfectly 
plastic  material  used  in  this  investigation.  These  can  be  stated  as  follows: 

(1)  It  is  assumed  that  the  Mises-Henchy  yield  condition 
accurately  determines  the  beginning  of  yield.  General 
considerations,  of  isotropy  and  symmetry  can  furnish  only 
the  general  form  of  the  yield  condition.  Beyond  this, 
any  yield  condition  is  a  hypothesis  which  only  tests  can 
justify. 

(2)  It  is  assumed  that  there  is  no  permanent  volume  change.  This 
assumption,  justified  on  the  basis  of  experimental  evidence 
for  metals,  leads  to  the  result  that  the  plastic  strain  is 
equal  to  the  plastic  deviator  strain. 

(3)  During  plastic  flow,  it  is  assumed  that  the  deviator  strain 
rate  tensor  is  proportional  to  the  instantaneous  deviator 
stress,  This  is  the  familiar  Prandtl-Reuss  postulate. 
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In  addition  to  these  three  main  assumptions,  it  is  possible  to  list  several 
other  restrictions  on  the  theory: 

(h)  The  material  must  be  isotropic.  This  condition  is  used 
in  developing  the  general  form  of  the  yield  condition, 

(5)  There  is  no  work  hardening,,  and  the  material  follows  the 
stress -^strain  diagram  of  Fig,  5  when  subjected  to  simple 
tension  or  compression. 

(6)  No  unloading  occurs,  Once  a  stress  point  has  yielded,  it 
remains  yielded  under  successive  increments  of  external 
load.. 

(7)  Time  effects  of  loading,  such  as  creep,  are  ignored. 

(8)  Displacements  are  small  so  that  the  small  deformation 
theory  of  elasticity  applies . 


3.2  Definitions  and  Notation 

The  following  definitions  and  notation  are  introduced  for  the 
purpose  of  describing  the  pertinent  constitutive  relationships  used  in  this 
work. 


T 

Total  Stress  Tensor  =  S  = 


a  T  T 

X  xy  ■ xz 

TO  T 

xy  y  yz 

X  X  a 

xz  yz  z 


(M) 


Spherical  Stress  Tensor 


s  0  Q 

0  s  0 


(15) 
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0 


s 
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Deviator  Stress  Tensor  =  S^  == 


S  T  T 

X  xy  xz 

T  c!  T 

xy  V 

T  ;  T  S 

xz  yz  z 


(16) 


where 


s  f=  mean  normal  stress  =  —  (a  +  a  +  a  ) 

3  ^  X  y  z' 


(IT) 


s  =  normal  x-component  of  S  =  a  -  s 

s  =  normal  y-comppnent aOf  =  a  ^  ^ 

y  y 

s  =  normal  z-component  of  =  (?  -  s 

z  z 


(18) 


With  these  notations,  it  is  obvious  that 


s  +  s  +  s  =0  +  0  +  0  -3s=0 

X  y  z  X  y  z 


Principal  nomal  stresses  are  designated  by  a^, 

Principal  normal  components  of  the  stress  deviator  are 
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(20) 
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Si  =  Si  -  s 
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S3  =  03  -  s 


A  completely  similar  notation  exists  for  strains . 
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Total  Strain  Tensor 
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e  0  0 


Spherical  Strain  Tensor  =  S'" 


0  e  0 


0  0  e 


Deviator  Strain  Tensor  =  E 


g  L  y  i  ^  ' 

X  2  'xy  2  'xz 

1  7  e  —  7 

2  xy  y  2  yz 


where 


1  7  —7 

2  xz  2  yz 


e  =  mean  normal  strain  =  —  (e  +6  +  e  ) 

p  X  y  z 


e  =  normal  x- component  of  E  =  e  -  e 

X  X 

D 

e  =  normal  y-cpmponent  of  E  =  e  -  e 

Y  y 

e  =  normal  z -component  of  E^  =  e  '  -  e 
'Z  z 


With  these  notations  it  is  obvious  that 


e'^  =  E®  +  E^ 


e  +e  +e  =€  +e  +e  -3e=0 
X  y  z  X  y  z 


Principal  normal  strains  are  designated  as  0^. 


Principal  normal  components  of  the  strain  deviator  are. 


ei  =  ei  -  e 


02  =  ^2  -  e 


e^  =  -  e 


(31) 
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3 • 3  Elastic  Stye ss -Strain  Relations 

In  the  elastic  range  the  relationship  between  the  elements  of  the 
stress  and  strain  tensors  is  assumed  to  be  that  of  Hooke's  law,  It  is  con¬ 
venient  to  express  this  linear  relationship  in  terms  of  the  elements  of  the 
deviator  stress  and  deviator  strain  tensors  as  follows; 


s 

=  2Ge 

S 

=  2Ge 

s 

=  2Ge 

X 

X 

y 

y 
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z 
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=  G7 

T  =  G7 
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xy 
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yz 

a  +  a  +  a  =  3K(€  +  €  +  e  ) 
X  y  z-^'x  y  z 


where 


G  = 


E 


2ll^^ 


^  3(l-2v) 

Eg.s . ,  '('32)(  can  bei  expressedimore  -concisely ■•as 


S°  =  2GE® 


(32) 


(33) 


(35) 


(36) 


Note  that  Eqs.  (32)  or  Eqs.  (36)  are  not  six  independent  relations;  since 

addition  of  s  +  s  +  s  =0  gives  an  identity.  Accordingly,  Eq.  (33)  is 
X  y  z  ' 

needed  to  give  a  complete  statement  of  Hooke ' s  law . 


3.k  Yield  Criterion 

A  yield  criterion  can  be  defined  as  a  condition  defining  the  limit 
of  elasticity  under  any  possible  combination  of  stresses  (8) ,  The  following 
considerations  of  isotropy  and  symmetry  show  what  the  general  fom  of  the 
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yield  criterion  must  te.  The  Mises  criterion  is  then  presented  and  reduced 
for  the  plane  strain  condition. 

Since  Hooke's  law  is  presumed  to  he  valid  in  the  elastic  range,  the 
strain  at  the  very  first  instant  of  plastic  deformation  is  \iniq_uely  determined 
by  the  stresses.  Thus  for  this  very  first  occurence  of  plastic  straining, 
the  yield  condition  can  be  written  as  a  function  of  the  stresses  alone. 


f(a  .  o,  T 


xy’ 


yz 


(57) 


Since  the  material  is  assumed  to  be  isotropic,  the  value  of  f  must  not  change 
if  the  coordinate  axes  are  rotated.  In  other  words,  f  must  be  an  invariant 
of  the  stress  tensor.  The  form  of  f  can  be  further  restricted  by  noting  that 
mere  hydrostatic  pressure  does  not  produce  appreciable  plastic  deformation  in 
metals  (8) .  Therefore  f  is  restricted  to  be  an  invariant  of  the  deviator 
stress  tensor. 

Let  the  deviator  stress  tensor  be  referred  to  its  principal  axes. 
The  following  three  linearly  independent  stress  invariants  are  then  chosen. 


Si  +  S2  + 


^  1/22  2,.  1/22  2,.  2  2  2 

=  —  (s^+s^+s,;  =  -  (s  +S  +S  j  +  T  +  T  +  T 
2  212  5  2  xyz  xy  xz  yz 


(58) 


J,  =  -v  (sf  +  S^  +  S'?) 
5  5  1  2  5 


Now.,  any  invariant  of  the  deviator  stress  tensor  can  be  expressed  in  terms 
of  these  three  linearly  independent  stress  invariants  (17) •  But  f  is  an 
invariant  of  the  deviator  stress  tensor.  Therefore  it  must  be  possible  to 
express  f  from  Eq.  (57)  in  terms  of  Jg  and 


J5) 


0 


(59) 
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The  yield  criterion  which  is  used  in  this  investigation  is  that 
of  Mises-Hencky  and  follows  the  general  form  above: 

Jg  -  =  0  (i^O) 

where  k  is  the  yield  limit  in  simple  shear.  Note  that  this  criterion  depends 
only  on  Jg.  For  equivoluminal  plane  strain  conditions,  Eq;.  (40)  reduces  to 

(-^)  -+  -  k^  =  0  (U) 

■  Eq,,.'  ''  ('%>)■  i(s:-'itiiebfdi5n-  df>:ithe>yite!Ldyc’C)ri!dliti'oh'  actualOiy.ius'e'd'  in;  the  .model  . 


3.5  Plastic  Stress-Strain  Relations 

In  order  to  relate  stress  and  strain  in  a  material  which  is  sub¬ 
mitted  to  plastic  flow,  it  is  convenient  to  express  the  strain  tensor  in 
terms  of  elastic  and  plastic  components.  Single  primes  will  be  used  to  denote 
an  elastic  component,  and  double  primes  will  denote  a  plastic  component.  Dots 
will  denote  rate  of  change  with  respect  to  increment  of  external  load. 

The  essential  nature  of  the  relations  between  stress  and  strain 
during  plastic  flow  is  contained  in  assumptions  (2)  and  (3)  of  section  3.1. 

The  assumption  that  there  is  no  permanent  change  of  volume  is  expressed 
mathematically  by  Eq.  (42). 

e"  =  i  (e;;  +  e;  +  e^)  =  0  (42) 

This  implies  that  the  plastic  strain  deviation  is  identical  to  the  plastic 
strain,  or. 
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Assumption  (5)  ojf  section  3 -I  states  that  during  plastic  flow  the 
deviator  strain  rate  tensor  is  proportional  to  the  instantaneous  deyiator 
stress  tensor.  This  is  expressed  mathematically  hy  Eg..  (^4)  helow; 


2Ge" 

=  Xs 

2Ge" 

;=  \s 

2Ge" 

=  \s 

X 

X 

y 

y 

z 

CG7'" 

G7  " 

= 

G7  " 

=  \T. 

,xy 

xy 

xz 

XZ 

yz 

(if4) 


where  \  is  a  proportionality  factor.  Eqs.  (44)  are^in  ithfe  same  form  as 
the  elastic  stress-strain  relations  given  in  Eqs.(32). 

The  hasic  relationships  which  are  assumed  during  plastic  flow  have 
now  "been  presented.  It  is  now  necessary  to  apply  these  relations,  along  with 
the  yield  criterion  Eg.  (40)  and.  the  elastic  relations  of  Egs.  (32)  and  (33) 
in  order  to  develop  the  final  relationships  hetween  the  stress  rates 
(incremental  stresses),  strain  rates  (incremental  strains),  and  the  instan¬ 
taneous  stresses. 

The  plastic  strain  rates  have  "been  expressed,  in  terms  of  stresses 
hy  Egs.  (44).  Similarly,  the  elastic  strain  rates  are  expressed  in  terms  of 
stress  rates  hy  differentiating  Egs.  (32): 
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(45) 


Combining  the  elastic  and  plastic  strain  rates  gives  the  total  strain  rate. 
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Note  that  these  relations  apply  only  during  plastic  flow,  i.e.,  when 


Jg  =  k  and  Jg  =  0 


(^7) 


In  order  to  eliminate  the  proportionality  factor  \  from  Eqs.  {,h6) , 
it  is  convenient  to  introduce  the  notation 


W  =  se  +  se  +  se  +  ‘r  j  +1:  7  +t  7 

XX  y,  y  z  z  xy  xy  xz  xz  yz  yz 


(48) 


where  W  may  he  interpreted  as  the  rate  at  which  stresses  do  work  during  a 
change  of  shape  and  to  note  that 


=  s  s  +  s  s  +  s  s  +  2t  t  -I-  2t  t  +  2t  t  (49) 

2  XX  y  y  z  z  xy  xy  xz  xz,  yz  yz  ' 


By  multiplying  the  first  three  of  Eqs.  (46)  hy  s  ,  s  ,  s  and  the  last  three 

X  y  z 

of  Eqs.  (46)  by  2t  ,  2,T  ,  2t  ,  respectively,  and  adding,  there  results 

.xz  yz 


•  .  2.2.2 

2GW  =  ss  +\s  +SS  +\s  +  SS  +\s 

XX  X  yy  'y  zz  z 

2  2  2 
+  2x  t  -I-  2\T  +  2x  T  +  2kT  +  2t  T  +  27.x 

xy  xy  xy  xz  xz  xz  yz  yz  yz 

i  2  2  „2^2„2s 

=  +  \(s  +  s  +  s  +  2x  +  2x  +  2x  ) 

2  '  X  y  z  xy  xz  yz' 

=12+  2XJ2 


(50) 


But  during  plastic  flow. 


Jg  =  k  and  Jg  =  0 


(47) 


Hence, 


2GW  =  2Xk 


(51) 


and. 


\  = 


GW 


(52) 
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Sulastituting  this  value  of  X  into  Eqs.  ,  it  is  pospihle  to 
solve  for  the  deviator  stress  rates,  which  gives  the  following: 


s  =  2G(e 

X  '  X 

s  =  2G(e 

y  y 

s  =  2G(e 
z  z 


_¥ 

2k‘ 

¥ 


2  "x) 


2k 

_¥ 

2k‘ 


2  "y* 


2 


-i  =  G(7 - o  ) 

xy  ^  xy  j^2  xy' 


T  =  G(7 - 3-  T  ) 

xz  '  xz  ,  2  xz ' 
k 


T  =  G(7  -  t  ) 

yz  yz  ^2  yz' 


(55) 


To  obtain  the  total  stress  rates  it  is  necessary  to  add  the  deviator 
stress,  rates  from  Bqs.  (55)  “to  the  spherical  stress  rate  which  can  be  obtained 
by  differentiating ,Eq.  (55) 5 


s  =  3Ke  (54) 

Adding -Eqs.  (53)  and  (54)  results  in  the  total  stress  rates,  as  follows: 


a  =s  +s=2G(e  -  s  )+5Ke  f  =G(7  ,  -  „ J 

XX  '  X  2j^2  x'  xy  ^  xy  j^2  xy' 


0  =s  +s=2G(e 

y  y 


s  ).+3Ke  ■i  =G(7  -  -4  ) 

y  2k^  y  xz  \xz  ^2  xz' 


(55) 


0z=s^+s=2G(ez  -  ^  sJ+3Ke 


T  =G(7  -  4  ) 

yz  yz  ^2  yz' 


Eqs.  (55)  give  the  desired  relationships  between  the  stress  rates,  strain 
rates,  and  instantaneous  stresses. 


5.6  An  Incremental  Form  of  the  Plasticity  Relations  for  Application  to 
the  Model 

In  general,  the  application  of  the  plasticity  relations  to  the  model 
is  closely  associated  with  the  three  stages  of  material  behavior  presented  in 
sections  3.3  through  3'5-  The  applications  of  Hooke's  law  and  the  Mises-Hencby 
yield  criterion  to  the  model  are  straightforward,  since  the  strains  can  be 
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computed  from  the  displacements  by  relations  similar  to  Eas.  (l)  and  the 
stresses  (or  forces)  at  a  stress  point  can  be  computed  directly  from  the 
displacements  by  relations  similar  to  Eqs.  (7).  Accordingly;  the  discussion 
which  follows  is  concerned  with  the  development  of  an  Incremental  form  of  the 
Prahdtl-Reuss  relations  for  application  to  the  model. 

Eas.  (55)  are  first  reduced  to  an  incremental  form.  Note  that  for 
plane  problems  the  number  of  relations  is  reduced  from  six  to  three. 
Therefore, 


Aa 

=  As  +  A  s 

X 

X  :: 

Aa 

=  As  +  As 

y 

y 

At 

xy 

=  G(A7 
'  xy 

AW  V 


For, plane  strain  conditions,  Eas.  (55)  are  reduced  to 


(56) 


-  fv 


and  Ea.  (5^)  beccmes 


As  =  5KAe  =  K(A€  +  Ae  ) 
'  X  y' 


The  incremental  W  becomes 


AW  =  s  Ae  +  s  Ae  +  s  Ae  +  T  A7 

XX  yy  zz  xyxy 


But 


S  =  a  -  rr  Ca  +  CT  +  Cf  ) 
z  z  5  X  y  z. 


(57) 


(58) 


(5^) 


(60) 


Where,  for  plane  strain. 
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a  =  V  (a  +  a  ) 
z  '  X  y 


(61) 


and  during  plastic  flow 
V  =  1/2 


(62) 


Hence, 


s  +  a)-  —  (a  +  a  +  =  0 

z2'x  y  j  x  y  2^ 


(63) 


Thus, 


^  a  +a  a  -a 

3  =  ^  .  i  (a  +  a  +  -V^)  =  ^ 

X  X  5  X  y  .  2  ^  2 


s  =  a  - 


(a  +  a  + 


—  \  IS/  I 

y  y  3  y  x  .2 


0  +0  a  -a 


-  -s 


(64) 


,  26  -e 

e  =  e  -  i  (e  +  e  )  = 

X  X  3  X  y'  5 
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Ae  i 


2€  -£ 

1  /  \  y  X 

e  =e  -  •^(e  +  6)  =  ■ — - 

y  y  3  X  y'  3 


2A€  -  Ae 

Ae  =— i- - ^ 

y  3 


Substituting  these  values,  Eqs.(64),  into  Eq.  (59)  yields  an  incremental  AW, 
which  reads. 


AW  = 


(a  - 

'  X 


a  )  (Ae 
y  '  X 


Ae  ) 

y 


+  T  A7 
xy'  xy 


(65) 


Using  the  expressions  for  AW,  and  from  Eqs,  (64)  and  (65)  in  Eq.  (56), 

Aa  becomes 

X 
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r  2A€  -  Ae 


Aa^  =  2G 


X  y 
5 


1  (a  -  a  ){Ae  -  Ae  )  +  t  A7  or  -  g  • 

2  ^  X  X.  y-  .xy  xy  ^  x 

2k^  ^ 


+  K(A€^  +  Agy) 


(66) 


Collecting  terms  gives 


Act  =  Ae 

X  X 


to  +  3K  _G_ 

^  k  • 


+  Ae 
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xy 
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(67) 


Similar  expressions  are  obtained  for  Act  and  At  .  as  follows: 

•7  ^*v 
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(69) 


These  last  three  equations  are  the  incremental  relationships  with 


which  the  incremental  stress  components  in  a  plastic  region  are  computed. 
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These  incremental  stress  components  are  added  to  the  existing  stress  components 

at  a  stress  point  to  obtain  the  total  stresses  acting  at  a.  yielded  stress  point. 

In  order  to  compute  the  quantities  ‘^^xy  appear  in 

Eqs.  (67)  through  (69),.  two  sets  of  displacements  corresponding  to  two  conr- 

secutive  load  levels  are  required.  One  set  of  displacements  is  the  set  which 

is  being  generated  for  the  current  level  of  external  load;  the  other  set  is 

that  computed  for  the  previous  external  load  level.  The  quantities  Ag  ,  Ae  , 

X  y 

Ay^  are  computed  as  the  differences  in  strains  determined  froin  these  two  sets 
of  displacements . 


IV.  SYSTEMATIC  RELAXATION  PROCEDURE  EOR  DETERMINING  DISPLACEMENTS 


4.1  Preljjginary  Remarks 

When  a  prohlem  in  continuum  mechanics  is  replaced  by  a  correspond¬ 
ing  problem  in  particle  mechanics  involving  a  discrete  model,,  the  question 
of  how  to  determine:  the  equilibrium  displacements  in  the  model  arises. 

Perhaps  the  most  obvious  solution  is  to  write  and  solve  the  set  of  simultaneous 
linear  algebraic  equations  (equations  similar  to  Eq.  (8))  for  the  unknown 
displacement  components  u^  and  v^.  Such  an  approach  has,  significant  dis.. 
advantages,  however.  The  preparation  of  the  equations,  whether  it  be  done  by 
hand  or  by  an  intricate  program  for  the  computer,  involves  a  considerable 
amount  of  labor.  In  addition,  even  with  machines  as  large  as  the  IM  709O 
the  number  of  equations  which  can  be.  solved  by  the  standard  library  subroutines 
is  limited  to  about  I50.  And  perhaps  most  important,  the  changes  in  the 
coefficients  for  the  displacements  resulting  from  yielding  of  one  or  more 
stress  points  are  not  at  all  easy  to  determine. 

A  more  flexible  and  practical  approach  to  the  problem  is  the 
relaxation  procedure  described  below.  Such  an  approach  eliminates  completely 
the  preparation  of  simultaneous  equations,  and  can  handle  a  very  large 
number  of  displacement  components  (of  the  order  of  several  thousand).  An 
additional  advantage  of  the  relaxation  method  is  the  physical  meaning  that 
can  be  attached  to  each  step  of  the  procedure.  This  is  of  considerable  help 
in  determining  plastic  stresses  and  strains. 

4.2  The  Relaxation  Procedure 

The  relaxation  procedure  used  for  determining  the  displacements 
can  be  graphically  summarized  by  means  of  the  flow  diagram  presented  in  Fig,.  6'. 
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All  mass  points  of  the  model  are  initially  in  equilibrium  with  zero  displace-r 
ments  and  no  external  load.  The  first  increment  of  external  load  is  then 
applied  to  the  boundary  mass  points  (or  other  specified  mass  points),  thus 
destroying  the  equilibri\jm  of  the  loaded  mass  points.  The  following  opera¬ 
tions  are  then  performed  for  each  mass  point  of  the  model. 

The  forces  acting  on  a  mass  point  are  determined  as  follows . 
■External  forces  acting  on  the  mass  point  are  given  as  a  part  of  the  loading 
pattern  applied  to  the  model.  Internal  forces,  originating  at  the  stress 
points,  are  determined  uniquely  in  the  elastic  range  from  the  displacements 
surrounding  the  stress  points  by  equations  similar  to  Eqs.  (7).  After  a 
stress  point  has  yielded,  the  force  components  at  that  stress  point  are 
determined  both  by  the  surrounding  displacements  and  the  past  history  of  thqt 
particular  stress  point.  Incremental  plastic  forces  are  determined  from  the 
incremental  plastic  stresses  given  by  Eqs.  (67),  (68),  and  (69).  These 
incfemental  plastic  forces  are  then  added  to  the  last  set  of  equilibrium 
forces  at  the  stress  point  to  obtain  the  current  total  plastic  forces  acting 
at  the  yielded  stress  point. 

After  the  forces  acting  on.  a  given  mass  point  are  determined,  a 
summation  of  all  the  forces  acting  in  the  x  direction  is  made.  In  general, 
this  will  result  in  a  residual  force  which  is  an  indication  of  the  amount 
by  which  the  mass  point  is  out  of  equilibrium  in  the  x  direction.  The  mass 
point  is  then  displaced  through  a  small  distance  in  the  x  direction  equal 
to  the  product  of  the  residual  and  a  flexibility  coefficient. 

Similar  operations  are  performed  for  the  y  direction.  This  places 
the  current  mass  point  in  equilibrium, . though  in  general  it  will  destroy  the 
equilibrium  of  surrounding  mass  points  by  a  small  amount.  The  procedure  is 
repeated  for  each  mass  point  until  every  mass  point  has  been  moved  once  in 
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the  X  direction  and  once  in  the  y  direction,  thus  completing  one  cycle  of 
relaxation . 

After  every  relaxation  cycle,  each  mass  point  is  inspected  to 
determine  if  it  is  in  eq.uilihritun.  If  not,  the  relaxation  process  is  repeated 
until  all  mass  points  are  in  equilibrium  to  within  a  prescribed  allowable 
error.  When  all  mass  points  are  in  equilibriim,  then  all  the  stress  points 
are  inspected  for  yielding  by  the  Mises-Henclqy  yield  criterion,  Eq.  (4l),,  and 
the  yielded  regions  are  recorded.  All  the  displacements  and  stresses  for  the 
equilibria  configuration  just  obtained  are  also  recorded.  If  desired,  the 
external  load  is  given  a  new  increment  and  the  complete  procedure  is  repeated 
for  each  load  increment  in  order  to  trace  the  development  of  plastic  yielding 
from  one  stress  point  to  another.  The  following  example  demonstrates  the 
manner  in  which  the  computations  are  perfo.rmed. 


4.. 3  A  Computational  Example 

Consider  the  elementary  example  shown  in  Fig.  7*  Only  mass  points 
"43"  and  "53"  are  free  to  move;  due'..to  symmetry  about  a  vertical  line  through 
these  mass  points,  the  u  and  v  displacements  at  a  mass  point  are  equal: 


"45 


"55  =  "53 


(70) 


Hence  there  are  only  two  unknown  displacements  in  the  problem,  Uj^^  and 

Using  the  material  constants,  dimensions,  and  loading  shown  in  Pig.  7^  it  is 

possible  to  write  two  simultaneous  linear  algebraic  equations  (similar  to 

Eq.  (8))  for  the  elastic  behavior  of  the  system  in  terms  of  the  two  unknowns, 

Ui ,  and  u^,.  Solution  of  these  two  equations  yields 

43  53 


'43  "  43 

U53  = 


-5  i. 

.252  X  10  inches 


(71) 
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Converting  .these  displacements  to  elastic  stress  components  by  means  of 
Eqs.  (7)  and  (4)  gives 

°x  "  ksi  a™  =  -.2143  ksi 

a®  =  -.0714  ksDi  =  -.0714  ksi  (72) 

y  y 

T  =-.2143  ksi  T  =  -.0714  ksi 

These  values  will  now  be  used  to  measure  the  progress  of  the  relaxation 
procedure . 

Before  beginning  the  systematic  relaxation  procedure^  it  is  first 
necessary  to  convert  external  pressures  to  concentrated  loads  for  application 
at  the  mass  points  and  to  determine  the  flexibility  coefficients  for  each 
mass  point.  For  example,  if  an  external  vertical  pressure  of  1  ksi  is  acting 
on  the  top  surface  of  the  model  shown  in  Fig:.  'J,  the  concentrated  vertical 
force  acting  on  mass  point  "43“,  which  arises  from  this  pressure  acting  over 
a  distance  of  \/2  =  I/2  inch  on  either  side  of  mass  point  "43”,  is 

=  (1  ksl)(|"  +  |")(1")  =  1  kip  (73) 

where  the  thickness  of  the  model  is  taken  to  be  one  inch.  This  vertical 
force  is  then  resolved  into  components  in  the  x  and  y. directions  for  applica¬ 
tion  at  mass  point  "43": 


P  =  .707  kip 

A 

Py  =  .707  kip 


(7^) 


The  flexibility  coefficient  for  a  mass  point  is  obtained  from  a  consideration 
of  the  effect  of  a  unit  force  acting  on  the  mass  point.  For  example,  a  unit 
external,  force  ■  of  one  kip  applied  in  the  x  direction  at  mass  point  "43"  is 
resisted  by  internal  force  components  acting  at  stress  points  "a"  and  ,"b." 
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■External  Load  =  1  kip  =  -(F^  -f  S 

X  xy 


(75) 


9.  1^ 

Expressing  F^  and  in  terms  of  displacements  by  means  of  equations  similar 
to  Eqs.  (7)  and  noting  that  all  displacement  components  except  are  beld 
fixed  gives. 


1  E  B  E  S 

^  "  (l+V')(l-2v}  S  ^  2  ■  2(l+v)  ^6^2 


(76) 


Solving  Eq.  (76)  for  u^^^  yields  the  flexibility  coefficient  in  the  x  direction: 


„  =  f  ^  lj-(l+v)(l-2v) 

X  (3--J+v)E 


(77) 


Because  of  the  symmetrical  arrangement  of  the  force  components  acting  on  mass 
point  "43;'j  the  flexibility  coefficient  in  the  y  direction  is  equal  to  f  • 


f  =  f 
y  X 


^•5  =  Ul+v)(l-2v) 
(3-4v)e' 


X 

(t8) 


A  similar  derivation  gives  the  flexibility  coefficients  for  mass  point  "53": 


f  55  _  ^  53  _  2(l+v)(l-2v) 
y  X  |.3-kv)E 

If  E  and  v  take  on  the  values  1000  ksl  and  0.25,  respectively, 
Fig.  7,  then  these  flexibility  coefficients  become 

f^^^  =  fy^^  =  .001250  inches/kip 

f^^^  =  =  .000625  inches/kip 


(75) 

as  shown  in 


(80) 


With  these  values  for  the  concentrated  external  loads  and  flexibility 
coefficients,  it  is  possible  to  begin  the  relaxation  procedure.  The  following 
step  numbers  make  reference  to  the  flow  diagram  of  Fig..  6. 
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Step 

1 

2 


3 

4 

5^ 

6 


7 


8 

9 

4 


Operation 

Set  u,  ,  =  V,  _  =  u^_  =;  v^,  =  0.  Also  set  force  components  =  0. 

■  45  43  53  53 

Apply  the  increment  of  external  load  to  mass  ppint  "43" • 

=  .707  kips 
Py  =  -707  kips 

Begin  with  mass  point  "43" ■ 

No  stress  point  has  yet  yielded,  since  aid  stress  components 
are  initially  =0.  Go  to  5h. 

On  the  first  cycle  all  force  components  are  computed  as 
zero,  since  no  mass  point  has  yet  teen  moved. 

On  the  first  cycle,  only  external  forces  ape  non-zero. 

Hence, 


=  +.707  kips 
=  +.707  kips 


New  =  old  +  f^''^ 

=  0  +  .00125  (-707)  =  .884  X  10"5  inches 
Similarly, 

=  0  +  .00125  (-707)  =  .884  X  10"^  inches 
Note  that  these  displacements  of  mass  point  "43"  destroy 
the  equilihrium  of  mass  point  ."53". 

The  current  mass  point,  "45",  is  not  the  last  mass  point. 

Go'  to  9. 

Take  mass  point  “55".  Got  to  4. 

Again  no  stress  point  has  yielded,  since  yielding  can  occur 
only  after  an  equilihrium  configuration  has  heen  reached. 


Go  to  5^. 


Operation 


Force  components  at  stress  points  "a"  and.  '’b"  are  computed 

from  Eqs.  (7)^  taking  account  of  the  eyanescence.  of  all 

displacement  conponents  except  u,  ,  =  v,  u^,  =  Note 

V  55  55 

that  only  those  components  acting  on  mass  point  "55''  are 
computed. 


y  ^  (l+y')(i-2v) 


(1-v)  -  V 


5 


1000 

(1+.25)(1-.5G) 


(l-.25)(o)  -  .25(.000881j.) 


=  -.177  kips 


V  ”  2(l+v) 
1000 

=  2(1+725) 


Im’ 

5  6 


0  -  .000884 


=  '•.177  kips 


X  '(l+v)(l-2v) 


(l-v)  ^  -V  ^ 


1000 

(i+.25)(i-.50) 


(l-.25)(0)  r  .25(. 000884) 


.177  kips 


E 


^xy  2(i+v) 
1000 


^  +  Z^' 

5  6 


2(1+.25) 


-  .000884  +  0 


=  -  .177  kips 


Note  that  for  the  first  cycle,,  mass  point, ''55,"  has  not  yet 


been  moved.  Hence  u^,  =  v^,  =  0  and  all  force  components 
55  55 
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Step 


7 


8 


Operation 

at  "m"  and  "n"  =  0.  From  considerations  of  symmetry,  it  can 
also  be  concluded  that 


S  =  S 
xy  xy 

The  equality  of  the  shearing  forces  and  the  axial  forces  at  a 
stress  point  on  this  first  cycle  is  purely  coincidental. 
Following  the  sign  convention  of  Fig.  2  for  positive 
forces , 


I 

I 


F  =  -S  -  f“  +  S  +  F^ 

X  xy  X  xy  X 

=  +  .177  +  "IT?  +  0  +  0  =  +  .554  kips 


F  =  -F®-S^  +  f’^  +  s'^ 
y  y  xy  y  xy 


=  +  .177  +  .177  +  0  +  0  =  +.55^  kips 


New  u^,  =  old  u^,  +  f 

53  53  X 


u^,  =  0  +  .000625  (-354)  =  .221  X  10"^  inches 
53 


Similarly, 

v^,  =  0  +  .000625  (.554)  =  .221  X  10"^  inches 
53 

Note  that  these  displacements  of  mass  point  "53"  destroy 
the  equilibrixom  of  mass  point  "43".  ” 

This  is  the  last  mass  point  and  the  end  of  the  first  cycle 
of  relaxation.  Go  to  10. 


10  All  mass  points  are  not  in,  equilibrium,  since  the  displace¬ 

ments  u.-,  and  v_,  under  step  7  above  destroyed  the  equilibrium 

53  53 

of  mass  point  “43" •  Hence  there  must  be  a  second  cycle  of 
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relaxation,  "beginning  at  step  3*  Note,  however,  that  in  only 
one  cycle  of  relaxation,  the  displacement  components  have 
attained  nearly  SO  percent  of  their  final  values. 

jThe  operations  listed  a"bove  demonstrate  the  procedure  for  elastic 
behavior.  Suppose  that  a  sufficient  number  of  relaxation  cycles  has  been 
performed  to  bring  both  mass  points  to  within  an  acceptable  error  in  the 
equilibriian  equations.  The  following  .discussion  indicates  how  the  yield 
criterion  is  applied  (step  11  of  Fig.  6)  and  how  the  force  components  at  a 
yielded  stress  point  (step  5a  of  Fig.  6)  are  computed. 

To  illustrate  the  application  of  the  yie.ld  criterion,  assume  that 
the  yield  stress  in  simple  tension  for  the  material  is  55  hsi.  Then  the  yield 
stress  in  simple  shear  is 

1,2  =  (-n|i^)  =  (^)  =  506  (-81) 


Applying  the  yield  criterion,  Eq.  (4l),  to  stress  point  "a”  gives 

_  306  <  0 

0.174  -  506  <  0 

and  to  stress  point  ''m“  gives 


(82) 


^-.214  +.071^' 


+  (-.071)  -  506  <  0 


0.010  -  506  <  0 


(85) 


Obviously  both  stress  points  are  far  from  yield  at  an  external  pressure  of 
only  1  hsi.  Indeed,  first  yielding  will  take  place  at  stress  point  "a"  at 
an  external  vertical  pressure  of 


-4o- 


=  k2  ksl' 


(84) 


Note  that  this  value  of  external  stress  is  considerahly  greater  than  the 
yield  stress  in  simple  tension  or  compression  of  35  ksi  ass\imed  for  the 
material.  This  is  characteristic  of  failure  or  yielding  in  two ' dimensional 
stress  systems,  and  will  he  evident  again  in  the  numerical  problems  presented 
in  Chapter  V. 

Until  the  load  level  has  reached  42,  ksi,  all  stresses  and  displace'- 
ments  increase  linearly.  When  this  elastic  limit  has  been  reached,  the  corre¬ 
sponding  displacements,  and  stresses  are  42  times  those  of  Eq.s.  (71)  and  (72); 

_2 

^45  =  =  4.242  X  10  inches 

u_,  =  V  =  1.058  X.  10  ^  inches 
55  55 

=  -55-00  ksi  =  -9.00  ksi 

X 

=  -5-00  ksi  =  -5. 00  ksi  (85) 

y  y 

=  -9-00  ksi  T  ™  =  -3,00  ksi 

xy  xy 


These  values  are  recorded,  and  are  used  to  determine  the  total  displacements 
and  stresses  for  the  first  load  increment  above  the  42  ksi  load  level. 

Suppose  now  that  the  load  level  is  increased  to  five  percent  above 
this  elastic  limit,  i.e.,  to  I.05  (42)  =  44.1  k&i,.  As  a  first  approximation 
to  the  final  displacements  at  this  new  load  level,  the  displacements  of 
Eqs.  (85)  are  also  increased  by  five  percent. 


^43  =  ’^43.  =  ^-^54  X  10 


-2 


U-,  =  v^,  =  1.111  X  10 

55  55 


-2 


inches 

inches 


(86) 


Note  that  two  sets  of  displacements  are  available :  the  last  set  of  equilibrium 
displacements,,  Eqs.  (85),  and  the  current  set  of  displacements,  Eqs.  (86) 


(which  in  general  are  not  compatihle  with  the  condition  of  equilibriirm) .  These 
two  sets  of  displacements  are  necessary. in  order  to  compute  the  incremental 
plastic  stress  components  according  to  the  discussion  in  section  5.6. 

In  order  to  compute  the  incremental  plastic  stress  components,  it 
is  necessary  to  compute  first  the  strains  at  the  stress  point  "a",  for  both 
levels  of  external  load,  by  Eqs.  (l) : 


•  For  load  level  =  h-2  ksi: 


a  .042^2 

X  8  1.U4 


-.05000 


a 

e 


y 


I53  ^  .01058 
8  1.414 


+.00749 


-y  ®  -  ’^33~^45  _  .01058  -  .04242 
^xy  S  1 . 4l4 

=  -.02252  (87) 


For  load  level  =  44.2  ksi: 


a 

e 


X 


.04454 

1.414 


-.05150 


a 

e 


y 


'..•01111 

1.414 


+.00786 


(88) 


.01111-. 04454 
1.414 


-.02564 


The  Incremental  strains  in.Eqs.  (67),  (68),  and  (69)  ere  obtained  by 
subtracting  Eqs.  (87)  from-Eqs.  (88): 


Ae^  =  -.05150  +.05000  =  -.00150 

Ae  =  .00786  -.00749  =  +.00057  (89) 

y 


-.02564  +.02252  =  -.00112 


it  is  convenient  to  compute  the  numerical 


Before  computing  Au  , 

X  y  xy 

values  for  G  and  K; 


E  1000 

°  -  2(l+v)  “  2(l+.25) 


400 


E  ^  1000 

3(l-2v)  3(1- -50) 


667 


(90) 


Note  that  instantaneous  values  of  the  stress  components  are  required  in 
Eqs.  (67);  (68),  and  (69)  in  order  to  compute  the  incremental  stress  com¬ 
ponents.  For  small  increments  in  the  external  loading,  the  instantaneous 
stresses  are  very  nearly  equal  to  the  stresses  at  the  last  equilibrium 
configuration, .Eqs.  (85). 

Substitution  of  Eqs.  (85),  (89),  and  (90)  intoEqs.  (67),  (68), 
and  (69)  gives  the  following; 


Aa^  =  (-.00150) 


4(400)+3(667)  400  ,- 


3  30^  ^  2 


(^) 


21 


+  (.00037) 


'-2(400)+3(667)  ,  400  /-33t3N 
3  30^  ^  2 


2-1 


+  (-.00112) 


-^0(-9)  p33+3n' 
'’”306  ^  2  ^ 


=  “.90  ksi 


(-.00150) 


-2(4oo)-f3(667)  4oo 
■3  ^  306 


(^)'' 


+ 


(.00037) 


4(400)-^3(667)  hOO  /-33+3')^ 
3  "  306  '•  2  ^  _ 


400 

30^ 


(91) 


+ 


(-.00112) 


45= 


=  (»..  00150) 


+  (.00057) 

+  (-.00112) 


.-4)0  /  ^\/-55+5\ 

•"5^  (-9)(^V^) 


H  (-9)(^) 


'<00  (1  -  4^) 


2  - 

w\ 


=  0 


Two  important  observations  can  be  made  immediately  from  inspection 
of  Eqs.  (91).  First,  the  stress  components  at  the  yielded  stress  point  "_a" 
are  not  increasing  linearly.  Second,  the  stresses  at  the  yielded  stress 
point  "a"  are  increasing  in  such  a  fashion  that  the  yield  condition,  Eq.  (4l), 
remains  satisfied.  This  is  a  consequence  of  the  fact  that  the  yield  condition 
is  used  to  eliminate  the  factor  of  proportionality  \  in  the  Prandtl-Reuss 
plastic  stress-strain  relations,  Eqs.  (44). 

To  obtain  a  first  approximation  to  the  stresses  and  forces  at 
stress  point  "a"  at  the  load  level  44.1  ksi,  it  is  necessary  to  add  the 
incremental  stresses,  Eqs.  (9l)l  to  the  last  set  of  stresses,  Eqs.  (85)? 

=  -.90-55.00  =  -55.90  ksi  4  =  4  I  "  ‘^5.95  kips 

a®  =  -.90-  5.00  =  -  5.90  ksi  =  o®  I  =  -  2.76  kips  (92) 

T  ^  =  0  -  9.00  =  -  9.00  ksi  s  ^  ^  =  -  6.56  kips 

Eqs.  (95)  correspond  to  step  in  Fig.  6,  wherein  the  forces  acting 
at  a  yielded  stress  point  are  computed.  Once  these  "plastic"  forces  are  known, 
the  relaxation  technique  proceeds  in  the  same  manner  as  before.  For  example, 
summing  forces  acting  on  mass  point  "45"  gives  the  result 

X^x  "  ^x'^^x'^^xy  ""  +31.20-25.95-6.56  =  +0.89  kip: 

)  F  =  P  +f\s  ^  =  +51.20-25.95-6.56  =  +0.89  kip 

Lj  7  yyxy 


(95) 


Hence  the  second  approximation  to-  the  displacement  of  mass  point  ."43"  is 
ol?tainod  by. adding  Eqs .  (86)  to  the  incremental  displacements  resulting 
from  the  unbalanced  forces  of  Eq.s.  (93)  ! 


Uj^j  =  .0445  +  .00125 (,89)  -  .0456  inch 

=  .0445  +  .00125(.89)  =  .0456  inch 


(94) 


where  .00125  is  the  flexibility  coefficient,  Eq..  (TT),  for  mass  point  "43''". 
Accordingly,  one  observes  that  the  displacements,  as  well  as  the  stresses, 
are  no  longer  linear  functions  of  the  external  load, after  plastic  yielding 


has  begun. 


V.  THE  NUMERICAL  PROBLEMS 


5 . 1  ProLlem  1;  A  Comparison  of  Theoretical  and , Model  Solutions 

ProLlem  1,  shown  diagraramatically  in  Fig.  8,  is  presented  in  order 
to  demonstrate  the  measure  of  accuracy  ohtainahle  with  the  model  used  in 
this  investigation.  The  theoretical  solution  is  obtained  from  that  given  hy 
Timoshenko  (2l)  for  a  single  concentrated  load  acting  vertically  on  the  sur¬ 
face  of  a  half-space.  To  obtain  the  approximate  theoretical  solution  for 
the  linearly  distributed  vertical  pressure  shown  in  Fig.  8,  the  effects  of 
seven  concentrated  loads,  located  symmetrically  with  respect  to  the  vertical 
center  line,  are  superposed. 

As  an  approximation  to  the  semi-infinite  half-space  of  the  theo¬ 
retical  solution,  the  following  boundary  conditions  are  used  for  the  model. 
The  left  boundary  is  assumed  to  have  a  zero  horizontal  displacement  and  a 
vertical  displacement  equal  to  that  of  the  material  spaced  a  horizontal 
distance  K  from  the  left  boundary.  The  lower  boundary  is  assumed  to  be 

I 

completely  fixed.  The  boundary  on  the  right  is  established  as  a  line  of 
symmetry.  These  boundary  conditions  are  indicated  graphically  in  Fig.  8. 

It  should  be  recognized  that  these  boundary  conditions  on  the  left  edge  and 
at  the  base  of  the  model  only  approximate  the  true  boundary  conditions  in 
the  half-space.  Accordingly,  exact  agreement  between  the  theoretical  and 
model  solutions  cannot  be  expected,  especially  in  the  regions  near  the 
boundaries • 

The  basic  solution  obtained  from  the  model  is  a  set  of  displace¬ 
ments  and  stresses  in  the  x  and  y  directions  oriented  as  shown  in  Fig.  1. 

For  presentation,  however,  all  displacements  and  stresses  are  resolved  into 
horizontal  and  vertical  components.  Figures  9;  10,  and  11  give  these 
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displacement  and  stress  components  for  Problem  1.  The  displacement  components 
within  a  square  refer  to  the  displacements  of  the  mass  point  located  at  the 
upper  left  corner  of  the  square.  The  stress  components  refer  to  the  stresses 
at  the  stress  point  located  in  the  center  of  the  square. 

To  facilitate  comparison  of  theoretical  and  model  solutions,  plots 
of  the  vertical  stresses  and  displacements  at  various  depths  in  the  half¬ 
space  are  given,  in  Figs.  12  and  IJ,  and  a  plot  of  vertical  deflections  at 
the  center  line  is  given  in  Fig.  l^i-.  Note  the  very  good,  agreement  of  the 
two  solutions  for  vertical  stresses  in  Fig.  12.  Only  near  the  lower  houndary 
is  there  any  observable  difference  between  model  and  theoryj  this  difference 
most  likely  arises  from  the  different  boundary  conditions  along  the  lower 
boundary  for  the  two  solutions.  The  pattern  of  vertical  displacements 
(Fig.  13)  appears  quite  reasonable,  and  the  comparison  of  these  deflections 
at  the  center  line  with  the  corresponding  theoretical  solution  (Fig.  l4) 
shows  a  good  agreement  in  the  pattern  of  the  deflections,  with  only  minor 
discrepancies  in  the  magnitudes  of  the  deflections.  Again,  this  difference 
in  the  magnitudes  of  the  deflections  obtained  from  the  model  and  from  the 
theory  of  elasticity  is  attributed  to  the  diffeFences  in  the  boundary  condi¬ 
tions  for  the  two  solutions,  particularly  the  condition  along  the  lower 
boundary . 

5.2  Problem  2;  Notched  Bar  Under  Tension 

As  an  example  of  a  type  of  problem  in  contained  plastic  flow  which 
can  be  solved  using  a  discrete  model  and  a  systematic  relaxation  procedure, 
a  bar  with  a  long  rectangular  notch,  or  slit,  is  shown  in  Fig.  I5.  In  the 
finite  model,  the  notch  actually  has  a  width  of  \,  though  for  practical 
purposes  the  notch  may  be  thought  of  as  having  infinitesimal  width.  A 


„4T» 

uniform  tension  is  applied  at  the  upper  edge  of  the  har,  the  left  edge  of 
the  bar  being  free  of  external  stress.  The  bar  is  assumed  symmetrical 
about  a  vertical  axis  through  its  center  and  symmetrical  about  a  horizontal 
axis  through  the  notch.  Hence  the  boundary  conditions,  on  the  right  and 
lower  edges  of  the  bar  are  those  of  zero  shear  on  the  boundaries  and  zero 
displacement  perpendicular  to  the  boundaries . 

As  mentioned  earlier,  the  basic  solution  obtained  from  the  model 
is  a  set  of  displacement  and  stress  components.  However,  once  successive 
sets  of  displacements  are  known,  the  stresses  can  be  computed.  Further, . it 
has  been  observed  that  the  general  pattern  of  stresses  does  not  vary,  appreci¬ 
ably  as  the  level  of  external  loading  is  increased,  even  though  portions  of 
the  material  may  be  undergoing  plastic  flow.  Accordingly,  only  the  basic 
solutions  in  terms  of  displacement  components  (Figs .  I6-I9)  are  given  for 
each  load  level  above  the  load  level  which  Initiates  plastic  yielding.  For 
this  elastic  limit  load  level  ^  complete  set  of  stress  components 

is  given  in  Figs.  20  and.  21,  and  plots  of  the  vertical  stresses  and  vertical 
displacements  for  various  depths  at  this  load  level  are  given  in  Figs.  22 
and  23. 

In  the  discussion  of  problems  in  contained  plastic  flow,  a  very 
useful  concept  is  that  of  an  "equiva.lent  shear  stress",  defined  as  follows; 

_  a  -Q  2  ^ 

Equiva.lent  Shear  Stress  =  j/jg  =  y  (--g  (95) 

Note  that  this  is  actually  the  largest  shear  stress  existing  on  any  plane 
passing  through  .a  given  point  at  which  and  are  defined.  If  this 

equivalent  shear  stress  is  divided  by  the  yield  stress  in  simple  shear,  k, 
the  ratio  represents  the  percentage,  of  the  yield  capacity  of  the  state  of 
stress  at  a  given  point.  Figures  2k,  25,  and  26  present  values  of  the 
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equivalent  shear  stress,  expressed  as  a  percentage  of  its  maximum  value  k, 
for  three  levels  of  external  load:  and  1.58a^^, 

It  is  of  some  interest  to  trace  the  development  of  the  yielded 
region,  as  the  level  of  external  load  Increases.  The  first  stress  point  to 
yield  is  the  one  at  the  very  end  of  the  notch  (Fig.  24).  It  is  of  signifi¬ 
cance  (Fig.  20  or  22)  that  the  vertical  stress  component  at  this  stress 
point  when  yielding  begins  is  47.4  ksi  --  considerably  greater  than  the 
assumed  yield  limit  of  55  ksi  in  simple  tension  or  compression.  As .mentioned 
previously,  this  is  characteristic  of  yielding  in  two-dimensional  stress 
systems;  the  yield  condition  depends  upon  a  combination  of  the  stress 
components  rather  than  on  the  value  of  any  single  component. 

To  be  strictly  correct,  the  external  load  increments  after  this 
first  stress  point  has  yielded  should  be  applied  in  very,  small  increments. 
Initial  investigations  indicate,,  however,  that  the  displacements  and  stresses 
are  very  nearly  linear  between  yielding  of  two  successive  stress  points, 
particularly  if  the  yielded  region  is  of  sm,all  extent.  Hence  the  next  two 
stress  points  were  yielded  by  relatively  large  increments  of  external,  load. 

At  an  external  load,  level  of  1.22a^^,  the  second  stress  point, 
immediate.ly  above  the  first  yielded  stress  point,  begins  to  yield.  As  the 
■load  is  increased  to  1.46a^^,  a  third  stress  point  yields  (Fig.  25).  Note 
that  the  yielding  is  not  taking  place  along  a  horizontal  line  at  the  waist 
of  the  specimen,  as  one  might  at  first  be  .led  to  expect,  but  is  progressing 
vertically  upward  and  to  the  right.  The  material  has  now  been  highly  enough 
stressed  so  that  only,  a  small  increase  in  externa,!  load  is  necessary  to 
propagate  the  yielded  region  completely,  across  the  bar  (Fig,.  26) .  In  problems 
of  this  type  which  involve  local  concentrations  of  stress,  the  specimen'  can 
actually  withstand  a  considerably  greater  external  stress  than  that  causing 
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initial  local  yielding.  Figure  27  summarizes  the  progression  of  plastic 
yielding  at  several  levels  of  external  load. 

This  pattern  of  plastic  yielding  shows  remarkably  good  agreement 
with  results  presented  by  Jacobs  (ll),  who  used  a  modified  stress  function 
approach  and  a  relaxation  technique  developed  by.  Allen  and  Southwell  (l) . 

As  Allen  and  Southwell  (l)  have  remarked,  this  type  of  plastic  yielding  may 
indicate  the  mechanical  behavior  behind  the  type  of  fracture  commonly  known 
as  "cup  and  cone".'  The  first  stages  of  failure  may  Involve  slipping  along 
planes  at  roughly  45  degrees  to  the  vertical.  Eventually  the  tensile  stress 
across  the  elastic  portion  of  the  waist  of  the  specimen  becomes  great  enough 
to  cause  a  breakdown  in  cohesion,  resulting  in  a  horizontal  tensile  fracture 
across  the  reduced  waist  of  the  specimen. 

Figure  28  illustrates  graphically  that  displacements  are  no  longer 
linear  functions  of  the  applied  loading  after  plastic  yielding  has  begun. 
■Load  deflection  ciorves  are  given  for  mass  points  located  at  "a",,  "b",  and 
"c"  of  Fig.  15.  .Mass  point  "a”  is  immediately .above  the  end  of  the  notch; 
mass  points  "b"  and  "c"  are  at  a  horizontal  distance  k/2  from  the  vertical 
center  line  and  at  vertical  distances  5"l/2  a>^d  2-l/2  k  from  the  horizontal 
center  line,  respectively.  Note  that  the  load  deflection  curves  differ, 
depending  on  the  location  of  the  mass  point,  and  that  the  load  deflection 
curve  for  the  material  within  the  elastic  core  at  the  center  of  the  specimen 
(mass  point  "c")  remains  nearly .elastic. 

5 . 5  Problem  5s  A  .Rirtially  Loaded  Half -Space 

As  a  second  example  of  a  problem  in  contained  plastic  flow,  the 
problem  of  a  partially  loaded  half-space  is  shown  in, Fig.  29.  Such  a  problem 
might  represent  the  effect  of  a  footing  on  soi.l,  or  a  machine  part  bearing 
against  another  part  of  much  larger  dimensions. 
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The  boundary  conditions  for  the  problem  are  the  same  as  those  for 
Problem.!;  and  the  elastic  solutions;  Figs.  50  and  through  is  quite 
similar  to  the  elastic  solution  of  Problem  1,  Preliminary  investigation  of 
plastic  yielding  under  the  triangular  loading  of  Problem,  ,1  indicates  quite 
■different  yield  patterns  for  the  two  problems;  however.  It  might  be  mentioned 
at  this  point  that  the  loading  pattern  shown  in  Fig,  29  purposely  introduces 
the  linearly  varying  stress  distribution  at  the  left  edge  of  the  loading 
pattern.  This  type  of  external  stress  distribution  reduces  significantly 
the  oscillation  in  displacements  and  stresses,  which  occurs  in  the  model  solu¬ 
tion  if  the  external  stress  distribution  drops  abruptly  from,  a  finite  value 
to  zero. 

The  concept  of  an  equivalent  shear  stress  is  again  used  as  a 
measure  of  the  closeness  to  yield.  Figure  38  shows  values  of  this  equivalent 
shear  stress  as  a  percentage  of  its  maximum  value  k  for  the  elastic  load 
limit  which,  initiates  plastic  yielding.  In  marked  contrast  to  the 

large  increments  of  external  lead  demanded  by  Problem  2  in  order  to  yield  a 
second  and  third  stress  point;  it  was  found  that  only,  a  small  increase  of 
two  percent  of  the  elastic  limit  load  was  required  to  initiate  yielding  at 
several  other  stress  points.  An  Increase  of  six  percent  (Fig,.  39)  in  the 
external  loading extended  the  yielded  zone  over  a  circular  arc  which 
almost  intersected  the  surface  of  the  half-space.  Figures  30"53  give  the 
basic  solutions  in  terms  of  displacements  for  each  load  level;  and  Fig,.  4l 
summarizes  the  progression  of  plastic  yielding  at  these  load  levels .  This 
pattern  of  plastic  yielding  under  a  partial  load  agrees  very  we.ll  with  the 
trajectories  of  maximum  shear  under  a  footing  given  by  Jurgenson  (12) . 

Note  again  (Fig.  36)  that  there  are  regions  within  the  material 
where  a  single  component  of  stress  (vertical  stress  immediately  beneath,  the 
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load;  for  example)  can  have  a  value  considerably  greater  than  the  yield  stress 
of  35  ksi  in  simple  tension  or  compression. 

The  non-linear  relation  of  load  and  displacement  at  specific  points 
within  the  material  is  also  evident  in  this  problem.  The  load-deflection 
curves  for  the  three  mass  points  "a";  “b",  and  "c"  of  Fig.  29  are  shown  in 
Fig.  ^l.  All  three  mass  points  are,  on  the  vertical  center  line;  "a"  is  at 
.the  surface;  and  "b"  and  "c"  are  at  depths  of  5A  and  8\  below  the  surface. 

The  surface  mass  point;  “a";  departs  greatly  from  the  linear  behavior;  since 
it  feels  the  curaulative  displacements  of  all  the  material  beneath.  Mass 
point  "b"  is  located,  within  the  yielded  zone  and  also  shows  a  non-linear 
behavior.  Maas  point  "e"  is  beneath  the  yielded  zone  and  exhibits  even  less 
than  linear  deflections .  This  seems  to  indicate  that  the  increments  in 
external  load  are  not  being,  transmitted  directly  through  the  yielded  zone; 
but  rather  are  being  carried  around  this  zone  by  a  redistribution  of  the 
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VI.  SUMMARY -AND  CON.CLUSI0N. 

The  object  of  the  thesis  is.  the  development  of  a  numerical  pro¬ 
cedure  for  the  solution  of  problems  in  contained  plastic  flow  of  plane 
continue.  To  accomplish  this,  a  discrete  model  is  introduced  to  replace 
the  physical  continuum.  The  equations  governing  the  behavior  of  the  model 
are  shown  to  be  identical  with  a  set  of  finite  difference  equations  for  the 
differential  equations  governing  the  plane  continuum. 

The  Mises-Hencty  yield  criterion  and  the  Prandtl-Reuss  stress- 
strain  relations  for  -plastic  straining  are  given,  and  a  finite  form  of  these 
relations  is  developed  for  application  to  the  model.  A  systematic  relaxation 
technique  for  the  computation  of  displacements  and  stresses,  within  the  model 
is  developed.  The  relaxation  technique  applies  to  both  elastic  and  plastic 
behavior,  and  is  well  adapted  for  use  on  large,  high-speed  computers. 

Three  numerical  example  problems  are  solved  by  means  of  the 
relaxation  procedure.  The  first  example  indicates  the  measure  of  accuracy 
obtainable  using  the  model.  The  last  two  examples  Illustrate  the  application 
of  the  procedure  to  problems  of  plastic  straining. 

Results  of  the  example  prcAjlems  indicate  that  the  numerical  pro¬ 
cedure  developed  herein  can  be  used  successfully  for -the  solution  of  a 
wide  range  of  interesting  and  practical  problems  in  contained  p.lastlc  flow. 
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FIG.  9  HORIZONm.  ARp  VERIIGAI  DISEEACEMENTS 
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FIG.  10  VERTICAL  STRESSES 
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FIG.  12  VERTICAL  STRESSES  AT  VARIOUS  DEPTHS 
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FIG.  16  HORIZONTAIi  MD  VERTIClAli  DISHACaiEKTS 
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FIG.  17  HORIZOUiaVL  AND  VEROUGAL,  DISHACHJ-IENTS 
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FIG.  18  HORIZONTAL  AND  VERTICAL  DISETACEMERTS 
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